A thermodynamic measure of the Magneto-electric coupling in the 3D topological insulator 
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We show that the magneto-electric coupling in 3D (strong) topological insulators is related to a second deriva- 
tive of the bulk magnetization. The formula we derive is the non-linear response analog of the Streda formula for 
Hall conductivity ( P. Streda, J. Phys. C: Solid State Physics, 15, 22 (1982) ), which relates the Hall conductivity 
to the derivative of the magnetization with respect to chemical potential. Our finding allows one to extract the 
magneto-electric coefficient by measuring the magnetization, while varying the chemical potential and one more 
perturbing field, a unique method never attempted before in the experimental search for the magneto-electric ef- 
fect. The relation we find also makes transparent the effect of disorder on the magneto-electric response, which 
occurs only through the density of states, and has no effect when the system is gapped. 

PACS numbers: 



The Quantum Hall effect 1 (QHE) was the first experimen- 
tal instance where a transport coefficient (the Hall conductiv- 
ity in d = 2) was quantized^. Finding an analog system in 
d = 3 had remained an unrealized dream of the condensed 
matter community, for many years. Recently, however, this 
has changed. The discovery of the topological insulatoi^l 
and specifically the 3D strong topological insulator (STTp^l 
have finally realized the dream of a d = 3 analog of the QHE. 
In the context of the QHE, Streda et alP^proved an insightful 
relation between the Hall conductivity, and a second deriva- 
tive of a thermodynamic potential 
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where p is the charge density, M the orbital magnetization 
per unit volume (perpendicular to the 2D system), B is the 
external magnetic field, and p the chemical potential. Here 
and throughout the manuscript we use units where c — e — 

h = \. 

Motivated by the analogy between the d = 2 QHE and 
d = 3 STI, in this paper we show that the magneto-electric 
coefficient characterizing the STI, is also related to a ther- 
modynamic derivative. We use an alternative formulation of 
magneto-electric effect, as a current response defined by 
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where cj> is an external field odd under inversion and time- 
reversal, and x is the magneto-electric response coefficient, 
and show that 
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There are a number of merits to this result, similar to those 
of ([TJ. First, |5]l suggests we can measure the Hall conduc- 
tivity in the QHE, by doing a thermodynamic measurement - 
vary the chemical potential through a back gate, and measure 
the magnetization of the sample. Calculating the derivative 
of the magnetization with respect to the gate voltage should 



give CHaii- Similarly, with ([3]) we can measure \ by simultane- 
ously varying the chemical potential p and the gradient (d a (f>), 
while measuring the magnetization of the sample. Second, 
the effects of disorder on ([T| and ([3) are entirely included in 
the density of states (DOS) D(e) through the particle density 
p = J def(e)D(e), where /(e) is the Fermi Dirac distribu- 
tion. Disorder will reduce the effective gap in the spectrum, 
compared with the clean limifiSI, but otherwise will not change 
anything for an insulator, as long as the chemical potential re- 
mains in the gap. 

The magneto-electric effecP^EH j s usually formulated as 
an anomalous term appearing in the action for the electro- 
magnetic fields in an insulator Seb — y~ J drdtP^E ■ B, 
with P3 the magneto-electric coefficient. Under inversion, E 
is odd, and under time reversal B is odd, so P3 is odd under 
both. Ref . fTTI showed that P3 takes on values modulo 1, and 
so in a material with either time reversal or inversion symme- 
try (or both), it can take on the values P3 = 0, h. The value 
P3 = \ then characterizes the STI. This value can in principal 
be measured^ as 
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where the Latin letters a,b — x,y, z denote spatial directions. 
Here and throughout the manuscript we will use the Einstein 
summation convention. 

At this point in time, a number of materials have been iden- 
tified as topological insulators using spectroscopy to charac- 
terize their unique surface state a 13 '' 17 l(an odd number of Dirac 
points). However, thus far measuring the magneto electric 
coefficient has proved challenging. The materials, by and 
large, have proven rather poor insulators, with significant car- 
rier concentration, in some cases even a bulk Fermi surface 
appearing 18 24 . In a bulk metal DC electric fields are screened, 
making it impossible to measure Q directly. 

Many other indirect ways have been proposed to detect 
the magneto-electric effect, most of them relying on the sur- 
face states in the ST j 11 | 12 | 25| BIl and on the Witten effec t 1 | 36 | 3 ^ l. 
The magneto-electric effect at the surface appears as a con- 
sequence of the material boundary, where P3 has a sharp 
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FIG. 1: Measuring the topological magneto-electric response with 
anti-ferromagnetic (AFM) order assuming the role of the auxiliary 
field cj>, odd under both inversion and time-reversal. At the micro- 
scopic level, some of the ions are non-magnetic (denoted by full cir- 
cles), while others are magnetic (denoted by open circles), and form 
AFM order (denoted by the arrows). A slow gradient in <f> in com- 
bination with the electromagnetic field generate a magneto-electric 
response. The <j> gradient indicated in the graph at the bottom is rep- 
resented in the image of the material by a change of shade (light red 
to light blue). 



jump, and the Witten effect appears as a consequence of a flux 
monopole entering the STL Our motivation in this work is to 
generalize the measurable effects of the magneto-electric cou- 
pling, beyond the surface and Witten effects, to response in the 
bulk. This can be achieved, as we explain below, by using the 
field <j) imitating the 4th vector potential component in 4D A4, 
and replacing P3. With this extra field <f>, we reformulate the 
magneto-electric effect with |2]), and find the thermodynamic 
relation |3]l similar to ([TJ. Formulated in this way, we will 
also be able to talk about magneto-electric effects in gapless 
systems, as was recently explored in Refs. 1381391 

Since the STI is the d = 3 analog of the QHE, we look 
to d = 2 for inspiration. First, the analog of magneto elec- 
tric coupling in a gaple ss metaP^Hs the intrinsic anomalous 
Hall effect (AHE) 40 " 43 ! in metals, where time-reversal is bro- 
ken, a large clean limit contribution to the Hall conductivity 
appears. The general expression for the intrinsic AHE con- 
tribution to the Hall conductivity becomes quantized once the 
chemical potential is in a gap of the band structure, giving 
the integer QHE. In a metal the effect of disorder cannot be 
neglected, and there are additional contributions to the Hall 
conductivity 43 . The interplay of the intrinsic AHE and disor- 
der can be understood using the so-called Bastin formula®4EI] 
for the electric conductivity tensor, which distinguishes be- 
tween contributions from states at the Fermi energy (the Fermi 
surface) a 1 and contributions from all filled states (the Fermi 
sea) a 11 . In an insulator, there are no bulk states at the Fermi 
energy, giving a 1 = 0. Streda et al. 18191 were able to show 
that the a 11 contribution can be related to a derivative of the 
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Here M, B are the orbital magnetization and magnetic field, 
respectively, ignoring Zeeman coupling to the electron spin. 
The formula holds in both gapped and gapless systems. For 
an insulator, a 1 = 0, and the total Hall conductivity is reduced 
to([TJ). 

The result ([TJ can be anticipated from the following con- 
siderations. In an insulator, with no dissipative currents, the 
only currents possible are persistent currents related to the or- 
bital magnetization J = V x M. The electric field is found 
from E = — Vjtt(x). Assuming the magnetization is an en- 
tirely local function of the intensive thermodynamic quantities 
M a = M a (T,B, M (x)),wefind 
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resulting in a xy — — ^Mf- . The second equality in |5]) is due 
to a Maxwell relatiorPEZl. 

As mentioned earlier, |5]l indicates the a 11 contribution can 
be found by doing a thermodynamic measurement. How- 
ever, the measured magnetization will include both the orbital 
and Zeeman contributions to the magnetization, while <|3j in- 
volves the orbital magnetization alone. For an insulator one 
can argue that the magnetization due to Zeeman coupling does 
not vary with chemical potential, and therefore measuring the 
derivative of the total magnetization, will give the same result 
as if we were measuring the orbital magnetization alone. For a 
metal on the other hand, the Zeeman effect magnetization can 
depend on the chemical potential, for instance in Pauli param- 
agentism. Therefore, measuring the total magnetization will 
only give a quantitatively accurate measure of a 11 in an insu- 
lator. Still, it will be useful in finding qualitative differences 
- it will exhibit quantization when the bulk is gapped. While 
such a measurement is conceptually straightforward, in prac- 
tice it is is more difficult than measuring Hall effect through 
electric currents. However, it has been carried ou i 48 l 49 l Next 
we will present the analog of |5]l in the d = 3 STI. 

In order to deal with a possibly gapless spectrum, as well as 
with disorder, we will have to formulate the magneto-electric 
coupling in a slightly different way from Refs. II 11121 The 
magneto-electric effect relates the magnetization to an applied 
electric field M = f^E equivalent to |4]). The Hall conductiv- 
ity is well-defined for both an insulator and a metal because 
it is defined through a current response to an external field 
Jx = PxyEy. For an insulator, J = V x M, and using Fara- 
day's law VxE = — <9 t B, we get 

J = ^VP 3 x E - ^8 t B . (7) 

Z7T Z7T 

The only DC (static) response comes from the first term 

J = ^VF 3 x E . (8) 
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Formulated in this way, we see P3 should be regarded as an 
external field. 

In fact, any inhomogeneous external field <fi with the same 
symmetry properties of P3, namely odd under inversion and 
time reversal, should suffice to generate such a current re- 
sponse. We therefore write J = ^^7<p x E, which is just 
Q. The response coefficient \ is the generalization of the 
magneto-electric coupling in the insulator, and a more faithful 
analog of the Hall conductivity - it can now be defined and 
calculated for metals as well 38 . For an insulator, \ w iU be 
quantized, yet the precise value of \ depends on how the field 
4> is defined and couples to the system, and therefore by it- 
self will not attain a universal quantized value. This is the one 
sacrifice we have to make in the alternate formulation of the 
magneto electric response. On the other hand, it will prove 
a more robust quantity to measure, in a system that may be 
gapless, and most importantly it will exhibit quantization in 
an insulator - the key qualitative feature we are after. 

Next we will derive ([3]). The derivation in the body of this 
manuscript is not rigorous, and does not apply to gapless sys- 
tems. It is presented here for the sake of brevity and clar- 
ity. In the supplementary material ?? we will derive the re- 
sult with some assumptions, while a general rigorous deriva- 
tion is left for a future publication. Using charge conservation 
dtp = —V ■ J, as well as Faraday's law VxE = — 9jB we 
find from Q 
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The second term from |7]) has no contribution since V • B = 0. 
Taking d c <fi = h c , the differential form, |9]) becomes 
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include any random potential. The chemical potential is in- 
cluded in H = ... — /i, so that the Fermi energy is at e = 0. 
The trace is over all degrees of freedom of the system - real (or 
momentum) space coordinates, and internal degrees of free- 
dom. Also, in the expression for the contribution x 11 > ai l 
Green's functions depend on the frequency e, while in the ex- 
pression for x 1 7 a U Green's functions have e = 0. Finally, 
the velocity is the conjugate operator to the auxiliary field 
H = Ho + f ^(x)?;^. The full details of this derivation 
we leave for a future publication. In the supplementary ma- 
terial ?? we provide a limited derivation, appropriate for an 
insulator, with the field 4> — h a x a coupled to a momentum 
independent v^. 

Much like in the case of conductivity, the form ( fTTj ) distin- 
guishes between Fermi surface contributions x 7 , which vanish 
for an insulator, and the contribution x 77 , which turns out to 
be the second derivative of the orbital magnetization, satisfy- 
ing <f3) 



where we used the same Maxwell relation as in d3). Note that 
the magneto-electric response is found from the magnetiza- 
tion parallel to the direction of the auxiliary field gradient. 
This should not be surprising as the magneto-electric effect 
should not care about whether the system is isotropic or not. 
In an isotropic system, the absence of any other directional- 
ity necessitates this outcome. Contracting the a, b indices in 
d = 3, we arrive at |3j. 

A rigorous proof of as well as a generalization to gap- 
less systems can be derived with some effort. Following stan- 
dard response theory techniques^, we can find the nonlinear 
response analog of the Bastin Formula 

Air 
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where v a are the velocity operators, /(e) is the Fermi Dirac 
distribution, S(x) is the Dirac delta function, and Gr = 
[e + iS — H] is the retarded Green's functions, which can 
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As we already discussed in the introduction, our result ( [12) 1, 
much like the Streda formula Q, indicates that x can t> e ob- 
tained by measuring the 2nd derivative of the (orbital) mag- 
netization. As noted earlier, at least in the insulating case, 
the Zeeman contribution to the magnetization should not vary 
with the chemical potential, and measuring the full magneti- 
zation instead of the orbital magnetization alone, will yield the 
same result. It is then conceptually straightforward to measure 
magnetization, and vary the chemical potential. The auxiliary 
field gradient h a , on the other hand, is at this point an ab- 
stract object we defined for our theoretical needs. We turn 
our attention now to discussing how <f> can be realized. First, 
given that <fi must be odd under time reversal and inversion, 
it can appear when anti-ferromagnetic (AFM) order is present 
in the material. It is not unimaginable that a topological in- 
sulator material could be stuffed with magnetic atoms that re- 
alize AFM order in the material. Second, we need <fi to vary 
(slowly) in space, as illustrated in Fig.[T] This can occur nat- 
urally in AFM order, as it tends to form magnetic domains. 
More difficult will be controlling and varying the strength of 
the AFM field. This can be done by changing the temperature 
of the system, and for better control of it, to be sufficiently 
close to the critical temperature of the AFM order. We men- 
tion in passing that if the material lacks inversion symmetry, a 
field breaking time reversal alone should suffice to generate cf>. 
However, a vast majority of the currently known topological 
insulators are inversion symmetric, and we therefore focus on 
this case in the current manuscript. 

We now turn to a concrete example for realizing as a 
spatially-varying Zeeman field in Bi2Se3 . The effective low 
energy continuum model derived for Bi2Se3 in Ref. |5U in- 
volves electrons in two orbitals. originating in the p-orbitals 
of different atoms (Bi and Se respectively). As a consequence, 
the two orbitals in general will have a different gyromagnetic 
ratio when coupled to a Zeeman field. Indeed, if the mag- 
netic field b is applied in the direction of the trigonal axis of 
the Bi2Se3 crystal (z-direction in the notation of Ref. [ST}, one 
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FIG. 2: Plot of x 11 /ffl numerically calculated for Bi2Se3 . The \ I! 
value is quantized as long as the chemical potential /i is in the bulk 
gap. Once /i is outside the gap, \ !I is no longer quantized. 

finds iceman = ba 3 (g + g 3 T 3 ), where 01,2,3 are the Pauli 
matrices of the electron spin, and t\ ^.3 are the Pauli matri- 
ces describing the orbital degree of freedom. If the magnetic 
field b varies on the length scale of atomic distances, it will 
effectively break inversion symmetry in the crystal and allow 
a more general Zeeman coupling to occur 

-^Zeeman — bFMO-3(go+93T3) + bAFM<^3(giTl+g2T2) ■ (13) 

Here bpM and bAFM are "ferromagnetic" (FM) and "antifer- 
romagnetic" (AFM) fields, respectively. The field bAFM is 
odd under both time-reversal and inversion, and is therefore 
a suitable realization of <fi. Generating it may require anti- 
ferromagnetic order, though any magnetic order that varies 
strongly on microscopic scales (ferrimagnetism, spin spirals 
etc.) will suffice. We also note in passing that we neglect the 
orbital coupling of the magnetic field we apply here. With a 
sufficiently weak Zeeman field the flux through a unit cell of 
the solid will be small, and we can safely neglect it. Using 
( fTTj i we calculate \ n for Bi2Se3 t with all the numerical pa- 
rameters we need (apart from gi) taken from Ref.[5T| Varying 
the chemical potential, we find the values plotted in Fig. [2] 
The most striking feature in plotting x 11 versus chemical po- 
tential is the plateau in its value while /i is in the gap. Once the 
chemical potential is outside the gap, the value of x 11 changes 
continuously. This will be the most easily discernible experi- 



mental signature of the topological insulator - a plateau in the 
transport coefficient x - even though its value is not universal. 

In conclusion, we have found that the magneto-electric cou- 
pling in topological insulators and their gapless counterparts, 
can be related to a third derivative of a thermodynamic po- 
tential. Most interestingly, this implies that the topological 
effects could be measured by probing either charge-density 
or magnetization in equilibrium, and in the bulk, rather than 
from non-equilibrium transport properties of the surface. Our 
result ([3]), suggests a conceptually simple way to measure 
the magneto-electric response, by measuring magnetization, 
while varying the chemical potential and the auxiliary field 
</>, or measuring the charge density while varying the mag- 
netic field and <j>. Moreover, our formula holds regardless 
of whether the system is gapless or gapped, clean or disor- 
dered. The measurement we propose, however, is challeng- 
ing. First and foremost, realizing the auxiliary field is diffi- 
cult, in the case of Bi2Se3 , requiring the introduction of mi- 
croscopic AFM order to the bulk of the material, and carefully 
controlling it. Controlling the chemical potential may also be 
challenging, given that we wish to probe 3D systems. Vary- 
ing the chemical potential is needed not only to calculate the 
derivative in Q, but also to detect the most clear cut evidence 
for a topological state - the plateau in x 11 > as illustrated in 
Fig. [2] Finally, the magnetization in our formula is the or- 
bital magnetization, ignoring the Zeeman contributions to the 
magnetization. In a metal, the Zeeman contribution can vary 
with the chemical potential, but in an insulator, it will not. 
Therefore measuring the full magnetization, instead of the or- 
bital magnetization, will yield x 11 in the insulating state, but 
in the metal it will yield x 11 P iu s some corrections. How- 
ever, the key qualitative feature is the plateau in x 11 m the 
insulating state, which will still show up when measuring the 
total magnetization instead of the orbital magnetization. It is 
the quantization of a response coefficient that signifies a topo- 
logical incompressible state. Despite these difficulties, our 
findings allow a unique conceptual approach to measuring the 
magneto-electric coupling, and it is our sincere hope this in- 
sight will be put to use in the lab. 
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